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Abstract 

We study the Ginzburg-Landau energy of superconductors with a term mod- 
elling the pinning of vortices by impurities in the limit of a large Ginzburg-Landau 
parameter k = 1/e. The function is oscillating between 1/2 and 1 with a scale 
which may tend to as k tends to infinity. 

Our aim is to understand that in the large k limit, stable configurations should 
correspond to vortices pinned at the minimum of and to derive the limiting ho- 
mogenized free-boundary problem which arises for the magnetic field in replacement 
of the London equation. The method and techniques that we use are inspired from 
those of [ pS3U (in which the case = 1 was treated) and based on energy estimates, 
convergence of measures and construction of approximate solutions. Because of the 
term as{x) in the equations, we also need homogenization theory to describe the fact 
that the impurities, hence the vortices, form a homogenized medium in the material. 

Key-words: Superconductivity, Ginzburg-Landau, pinning, homogenization. 
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I Introduction 



Superconducting materials have the property of expelhng an apphed magnetic field. In fact, 
the behaviour of a superconducting sample varies according to the value of the applied field 
and the value of the Ginzburg-Landau parameter k which is characteristic of the material. 
When K is large, the superconductors are known as type-II and display vortex patterns 
for intermediate fields: for high magnetic fields, the material is normal and the magnetic 
field penetrates into the sample, for low fields, the material is superconducting, that is the 
magnetic field is expelled from the sample and for intermediate fields, there are vortices. 
The vortex state is a state where the superconducting and the normal phases coexist: at the 
center of the vortex, the material is normal and the vortex is circled by a superconducting 
current carrying a quantized amount of magnetic fiux. The motion of vortices generates an 
electric field hence energy-dissipation. In order to have the desired property of dissipation- 
free current fiow, the vortices have to be held fixed or pinned. In practice, attempts are 
made to pin vortices either by varying the thickness of the material or by introducing 
impurities or normal inclusions. Sufficiently strong pinning is necessary for functional 
superconductors capable of sustaining strong currents and high magnetic fields. The new 
high-temperature (high Tc) superconductors are strongly type-II superconductors, that is 
their phenomenology is dominated by the presence and properties of vortices when an 
exterior magnetic field is applied. The pinning problem is particularly intricate in high-Tc 
superconductors where it depends on specific structures such as layering and structural 
defects. 

In this paper, we will be concerned with the case where the vortices are pinned by 
impurities in the framework of the Ginzburg-Landau model. We will study the behaviour 
of global minimizers of the Ginzburg-Landau energy when a term modelling the pinning 
of vortices by impurities is added, in the limit of a large Ginzburg-Landau parameter k, 
which describes extreme type-II materials. 

I.l The Ginzburg-Landau model with a pinning term 

Recall that in the framework of the Ginzburg-Landau theory (see fTj for more details), the 
state of the material is completely described by a vector potential A and a complex- valued 
function u, which can be thought of as a wave-function of the superconducting electrons, 
and is nondimensionalized such that \u\ < 1. The type of material is characterized by 
the Ginzburg-Landau parameter k and in the case of type II, k is large so that we define 
e = I/k, which will be small. The energy is the following: 

(LI) Mu,A) = ^ j \(y-iA)u\'' + ^^{a,{x)-\u\''f + \h-K^\\ 

n 

Here, Q is the domain occupied by the superconductor, h = curl A is the magnetic field and 
/lex is the exterior magnetic field which is constant in our problem. A common simplification 
is to restrict to a two-dimensional problem corresponding to an infinite cylindrical domain 
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of section i7 C (smooth and simply connected), for an applied field parallel to the axis 
of the cylinder. Then A : Q M."^, h is real- valued and all the quantities are translation- 
invariant. 

The energy that we are going to study here is slightly different from the classical 
Ginzburg-Landau energy in the sense that there is a term penalizing the variations of the 
order parameter u. We denote this function by ^^(x). In the case originally studied by 
Ginzburg and Landau, = 1. In this paper, a typical example for would be to oscillate 
between 1/2 and 1 in the domain, with a typical scale rj which may tend to with e. The 
minima of correspond to the impurities in the material. Hence it is expected that these 
minima will be the pinning sites for the vortices. 

The modified Ginzburg-Landau functional (|I. 1|) was first written down by Likharev ||L| . 
Then, this model has been used and developed in ||CK|| and ||CDG|| . Review articles on the 
topic include PFGLV |, |CT|, [ 02 ] and Q. Computational evidence that the vortices are 
attracted by the impurities, that is the points of minimum of a^ix) can be found in | |GDG | 
or [pUPj . 

In this paper, we want to address the question of how the term will modify the 
properties of the superconductor in the presence of an exterior magnetic field. The method 
and techniques that we are going to use are inspired from those of |pS3|] (in which the 
case Oe = 1 was treated) and based on energy estimates, convergence of measures and 
construction of approximate solutions. Because of the term a£{x) in the equations, which 
can be a rapidly oscillating function, we will also need homogenization theory ( ||CD|| , [|JKO|| , 
MuT|| ) to describe the fact that the impurities, hence the vortices, form a homogenized 



medium in the material. 



1.2 The equation for the magnetic field 

The Ginzburg-Landau equations associated to the functional ( p]) when minimizing for 
{{u,A) e H\nX) X H\n,M.'^)} are 

(GL) I ~ iAfu = ^uia^ix) - \u\^) 
\ -V^h =< iu, (V - iA)u >, 

with the boundary conditions 

h = /lex on dQ 

(Vm — iAu) ■ n = on dQ. 

Here V"*" denotes {—dx2, dxi), and < z,w >= Re{zW) for z, w in C. Recall that the problem 
is invariant under the gauge transformations 

u — s> ue 

where $ G if^(f2,M). Physically meaningful quantities are gauge invariant. These include 
the energy J^, the magnetic field h and the superconducting current j =< iu, {V — iA)u >. 
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Let us describe the properties of a superconductor. These phenomena are described 
for instance in Q. The state of the material depends on the apphed field /lex- In the 
absence of pinning, that is when = 1, there are two critical fields ifci and for which 
a phase transition occurs. Above ifc2 = 0{\), superconductivity is destroyed and the 
material is in the normal phase {u = Q,h = hex)- Below Hc^ = 0(|log e\), the material is 
superconducting everywhere, that is |m| ~ 1. This is the Meissner phase characterized by 
complete expulsion of the magnetic field : in the limit when e goes to zero, the magnetic 
field satisfies the London equation 



[L2) 



-Ah + h = inQ 
h = hpx on dQ. 



Between Hc^ and Hc2, the material is in the mixed phase defined by the coexistence of 
the normal and superconducting phases in the form of vortex filaments: the magnetic field 
penetrates into the material in the form of flux lines at the center of which u vanishes. The 
induced magnetic field approximately satisfies 



[L3) 



—Ah + h = 271 di6p. in Q 
h = hex on dfl, 



where the p^'s are the centers of the vortices, and the cij's their degrees, that is the topo- 
logical degree of the map u/\u\. These filaments are of characteristic size e. They are 
surrounded by a superconducting region in which |m| ~ 1. In order to minimize their 
repulsion, the flux lines form a triangular lattice, called the "Abrikosov lattice". With 
increasing fields, the density of flux lines increase until the vortices overlap and is 
reached. The generation of vortices by the external field has been mathematically studied 
very recently in 0, |S2|, ^, [SS^, [SS^. 



In ||SS3|| , it is proved among other things that, in the limit when e tends to 0, equation 
(|]3D is replaced by 



(1.4) — A/i* + /i* = /i* 

where fi^ is the density of vortices in units of hex and h^ = h/hex- The measure /i,,, is 
supported in an inner region u depending on the value of hex and is of uniform density in 



Our aim is to give a rigorous proof that in the small e limit, stable configurations should 
correspond to vortices pinned at the minimum of and to derive the limiting homogenized 
free-boundary problem which arises for the magnetic field in replacement of the London 
equation (|L4D . 

Using the second equation in (G.L.), we notice that the energy can be rewritten 

(1.5) Js{u,A) = ^ J ^^\Vh\'+\h-hex\' + ^ J |V|«|p + ^(a,(a;)-|np)2. 

n n 
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We will show that for a sequence of minimizers {u£,As), the second integral in ( (TSP is 
negligible. Then, when e tends to 0, j-up ~ ^^(x) outside the vortices, and our main result 
will state that /i^ = curl satisfies roughly the following equivalent of ( [TSD in the case of 
pinning: 



(1.6) -div (^^V/i,^ +K = 27^ J2d,Sp^. 



The existence of pinning will modify the locations Pi of the vortices and the value of i^ci- 
Since is a rapidly oscillating function describing impurities, the framework for passing 
to the limit when e is small is that of homogenization theory. When passing to the limit 
in (|1.6[), we obtain a different limiting operator from ([T^), that is 

(1.7) -div (AoVK) + K = fx^ 

where /x^, is a positive measure which is supported in an inner domain and .Aq is the 
homogenized limit of the matrix = — I in the sense of iJ-convergence, see definition 
below. 

Definition 1 We say that the family of 2 x 2 matrices A.£ H -converges to when e 
tends to 0, if and only if, for any f in H~^{fl), the solution in Hq{Q) of 

-diY{AeVVe) +Ve = f 

satisfies 

Ve^VQ weakly in Hq{Q), 
A^Vvs^AqWvo weakly in (L^(fi))^, 

where Vq is the Hq{Q) solution of 

-div{AoVvo) +vo = f. 

We refer to the work of Murat and Tartar ||MuT | for more details on the notion of H- 

convergence; one can also see ||CD| , [JKU|] . In the following, we will always let A.^ = — I- 

Then .Aq is also a diagonal matrix. In the general case, the computation of .Aq is hard 
and not always known, see | JKO | for examples. But in some simple cases, this definition 
allows to compute A^. For instance, if as{x) = a{x/e), and a{x) = ai{xi)a2{x2) where ai 
and 02 are periodic, then 



= diag(^,\\ with 



-/ 1 

a; = a,; I — 



where al denotes the mean of over a period (see ||JKO|| ) . Note that even though the 
sequence has no pointwise limit, the limiting problem and Aq are well defined. 

An important property of if-convergence (see |[Mu'l| ) is that if the sequence is 
bounded from below and above by positive constants independent of e, then there exists a 
subsequence A.^' and a matrix Aq for which As' if -converges to Aq. For us, it will imply 
in the following that up to the extraction of a subsequence, the family A^ ii-converges to 
some limit Aq, thus leading to the limiting problem (0). 
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1.3 Main results 

Let us now state our hypotheses and results. We assume that /lex is a function of e and 
that the following limit exists and is finite: 

(1.8) A = liml^°^"l 



£-0 h^^{e) 

Moreover, we make the following hypotheses on the function ae{x): 
(HI) There exists a constant 6o > such that 6o < dei^) ^ 1- 

(H2) There exist a constant C and a sequence ri{e) (which may tend to with e) such that 

C 

l/ri{e) < /lex and \Va^\ < ——. 

r]{e) 

(H3) There exist a continuous function b{x) and a nonnegative functions /^^(x) such that 
a^i^x) = h{x) + /3e(x) and for any e > and any x G fi, v[im.B[x,&{e)) Pe = 0, where 

5{e) < 



(log |log £|)2 

(H4) The family of matrices A-e i^-converges to >Ao- 

Note that, as we mentioned earlier, it follows from hypothesis (HI) and the compactness 
of the set of matrices bounded from above and below that there exists a subsequence of 
which i7-converges to ^4o ||MuT|| . Our hypothesis (H4) is there to restrict to this sub- 



sequence for ease of notation and to impose that the whole sequence converges. Moreover, 
(H2) means that can be a constant independent of e but can also oscillate very quickly 
with e (but not too quickly, i.e. not quicker than hex)- Note that in the case where does 
not depend on e, then A.^ = A« is constant. 

Let us emphasize that because /?£ > 0, 6 can be thought of as the lower envelope of 
and the local minima of are the local minima of h. Hence h will be related to the pinning 
sites of vortices and the oscillations of are those of /?£. Moreover, the hypotheses imply 
that b > bf). 

First, let us state the result concerning the limiting problem ([TTD. We relate h^, and fi^, 
to the minimum of a variational problem. Let Ai denote the space of Radon measures in 

n. 

Theorem 1 Let us assume that (HI) to (H4) are satisfied. Let us define for any A > 0, 
(1.9) Eif) = b{x) |-div(AV/) + f\ + ljvf- A)V/ +\f-l\\ 
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over 

V = {f s.t. / - 1 e Hlip), and - div(A)V/) + / G M}. 

The minimizer of E over V exists and is unique. It satisfies 
( K-le H^{Q) 



. . ^ 

Ai* > 1 — m \l 

/i* (k-{i-^)] =0 in n. 



Moreover yU* > and E H ^{^l). 

Problem (P) is a free-boundary problem, called in the literature an "obstacle problem" 
(see [0). Another way of considering problem (P) is to define the subset of Q 

(I.IO) uJA = {x e n, s.t. K = l- Kb/2}. 

Then /x^ = in f2 \ cJa", and h^^ = 1 — A6/2 in a;A, dujf^ being called the "free-boundary", 
because cja is unknown and uniquely determined by the set of equations (P). 

Note that if and h are smooth enough then h^, is C^'" (a < 1), /i* is in L°°, the 
free-boundary duj\ is regular for almost every A (see | |BlVi| | ) and then we can write 

, A6 A , . . 
/i^ = 1 - — + — div(>^V&) m wa. 

Once we have proved Theorem ^ concerning the limiting problem, we can get convergence 
for any sequence of minimizers (m^, A^) of the energy Je{ue, Air) to E{h.^) in a sense similar 
to F-convergence. 



Theorem 2 Let us assume that one? (^i^ij to (^ii/^y' are satisfied. Let ('Ue,^^) 6e a 

family of minimizers of J^, and hs = curlA^ the associated magnetic field. Then, as e 
tends to 0, 

— > weakly in H^{Q), 



cx 



where h^, is the minimizer of E. Moreover, 

(1.11) hm = E{h,) = ^ [ + \ !vh,- AoVh, + \K-l\\ 

(1.12) ^ > V/i* ■ ^V/i* + Afeu^K, m t/ie sense of measures. 

One can easily notice that if A = (i.e. if hex ^ |log e\), the solution of (P) is h^, = 1, and 
E{h^) = 0. In this case. Theorem ^ asserts that 

— > 1 strongly in H^, and lim — -^r— ^ = 0. 

hex e-*0 hi^ 

The proof of Theorem 2 is the main part of the paper (see Section 1.6 for a sketch). 
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1.4 The case A > 



Let us now present some stronger results in the case where A is positive, i.e. /lex is of the 
order of |log e\. The first issue is to determine mathematically the location of vortices. 
From the physics, we know that vortices are the zeroes of Us with non-zero winding number. 
Instead of defining vortices, we isolate them in disjoint vortex balls covering the set where 
is small. The centers of these balls can be thought of as being the centers of the 
vortices. 

Proposition I.l Let us assume that A > and that (HI) to (H4) are satisfied, then there 
exists Eo such that if e < Eq and {us, A^) is a minimizer of J^, there exists a family of balls 
of disjoint closures (depending on e) {Bi)i(zj^ = {B{pi,ri))i!zi^ satisfying : 

(1.13) [xen, \,/^-\u,{x)\\>-^\ C \ \B{p,,r,) 

(1.14) J]n< ^ 



1 f IV/? P 



|2 



Bi 



where = curlA^, and di = deg{j^^,dBi) if B^ C Q, and otherwise 



This proposition will be proved at the beginning of Section II. Here is the meaning of the 
different inequalities: ( |I.13| ) locates the set where Im^I differs from a^, which is contained 



in a union of disjoint balls; these balls represent the vortices or clusters of vortices. (P!.14| ) 



gives a control on the size of the balls and ( |1.15| ) gives a lower bound on the energy, which 



is the contribution of vortices according to their degree di and their location pi, appearing 
through the value b{pi). As opposed to the case of = 1 (see | fSS3| ), the least energy is 
attained for pi at the minimum of b. 

Using this proposition. Theorem can be made more precise: 

Theorem 3 Let us assume that A > and that (HI) to (H4) are satisfied. For any balls 
B{pi,ri) and integers di which satisfy ( \I. ISj )- ( fO^ j - then 



(1.16) lim — "^ diOeipi) = f^b\fx^\ 

(1.17) ^Y.d^Sp 

jG/e 

(1.18) l-El^^l'^^ 



in the sense of measures, where 
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1.5 Physical interpretations and consequences 

Our results show that /i*/iex is a good approximation of /i^ and that, in the hmit e 0, 
the vortices are scattered in an inner region uj\ with density fi^,, where = 1 — Ab{x)/2. 
In the outer region Q \ uJX, there are no vortices and satisfies —div{A^VK) + = 0. 
Unhke the case = 1, the vortex-density in uJX is non- uniform in general. Moreover, as 
A decreases, the vortex-region first appears at the minimum of ip as defined by problem 
( [1.19|) below: as in |^S3|| , we can derive a necessary and sufficient condition for to be 



nonempty. 

Proposition 1.2 Let ip be the solution of 



(1.19) 

then 



-div{Ao'Vilj) +11^ = -1 inn 

ijj = on Oil, 



hex 1 
u\ <^=^ lim — j- > 



e-^o I log e I 2 max 



If we define Hc^ as the field such that for h^x < Hci, the minimizer of the energy has no 
vortex (i.e. \u\ > &o/2) and for h^x > Hci, there exists a minimizer with vortices; then 



Proposition gives a hint that 

log el 



-'-'ci 



2 max \^\ 

Thus the presence of pinning modifies the values of the first critical field (see pSl 



for the case without pinning). In fact, we could adjust the proof of ||SS1|| to obtain: there 
exists ke = 0(1 log | log£:||) such that for e small enough and 

2 max Itpl ^ 

then any minimizer has no vortex. 

Furthermore, the position of the minimum of ip depends on the pinning potential a^ix). 
As A further decreases, corresponding to hex increasing, the vortex-region cja grows, until, 
for A = {hex ^ |log e\), cua = Q. At this point there are so many vortices that the 
macroscopic density of vortices and the induced magnetic field are no longer infiuenced by 
Og. In other words, the strength of fiux pinning is for hex ^ I loge|. 

In the case where a^ix) = a{x) is independent of e, a{x) = b{x) and ^4o = a^^T. Hence 
the limiting problem is a London equation with weight. We would like to point out that it 
is natural to define a vortex velocity by t> = j^2'^h (see ||(JyP| ] ) . In particular 



a 
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can be defined as a limiting velocity (per unit of hex)- Note that in uj\, since = 1 — ^Aa, 
then = — |AVloga. It implies that when a is constant, v^, = and there is no mean 
current in the vortex region. But when a varies spatially, there is a nonzero limiting mean 
current and a nonzero limiting velocity w*. Hence v ~ /lex'^* that is ^log/tVloga. This 
is the result of Chapman-Richardson ||CR|| in the case where the three-dimensional vortex 



line has no curvature. They describe the phenomenon saying that the variation in a acts 
as a pinning potential. 

When A = 0, the velocity v^, is zero as well. Decreasing A means increasing the field. So 
when a varies spatially, there is a critical exterior magnetic field above which the pinning 
potential has no role and the current is destroyed. 

In the general case where depends on e, it would be interesting to prove a convergence 
of the mean vortex velocity = ji^Vh^. Still, one can observe two different effects coming 

from the presence of pinning in the term \'Vh^\'^/a£ and resulting in the energy E{h.^) in 
the homogenization process: 

- One effect is related to the concentration of energy in the vortices and the location of 
the vortices. It appears through the term 

A 



2 



in the limiting energy E. This term is smaller if /i,. is non-zero at points where h is minimal. 
( [1.16 ) implies that vortices go to points where /3e = 0. These points will be called pinning 



sites in the following. Because 5{e) tends to 0, the number of such points is big. The effect 
on the position of vortices is to see h and the minima of h. Moreover, since ( [1.17|) and (P^.18| ) 
have the same limit, it means that vortices tend to have positive degrees. 

If h does not depend on x then /i* and /x* are constant in wa, and there is no change 
for the location of vortices from the case = 1. On the other hand, if h is non-uniform, 
then V/t* is non-constant in cua and there is a pinning current. If for example the domain 
is a disc and the minima of 6, that is the impurities, are located at sites different from the 
center of the disc, one expects that vortices, or the vortex-region cua will be closer to the 
minima of 6, but it seems difficult to give a rigorous proof of this qualitative fact. 

- The other effect is due to the rapid oscillations of with e and comes from the energy 
outside the vortices, converging to the homogenized term 

i/v/...Av..+i/.-ir^ 

in E. It changes the equation for the magnetic field h from the usual London equation. 
If /?£ 7^ 0, then the homogenization effect can be anisotropic. The size 5{e) (which can 
be related to t] if /J^ is not identically 0) cannot be taken bigger than in (H3), otherwise 
each pinning site would be too large and the vortices could push one another outside the 
pinning site. 
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Let us also point out that we cannot allow stronger oscillations of than in (H2), be- 
cause the second integral in ( jlT^ ) would become the dominant term. It would be interesting 
to investigate what happens if (H2)-(H3) are relaxed. 



1.6 Main steps of the proof 

Let us now state the two steps of the proof of Theorem ||. It is obtained as in PS3|| by 
getting first a lower bound on the energy, Proposition |L3| , proved in Section II, and then 
an upper bound. Proposition J^, proved in Section III. 

Proposition 1.3 Let us assume that A > and that (HI) to (H4) are satisfied. Let 
{us, A^) be a minimizer of J^. Then 

(L20) liminf ^ J^K, ^) > ^ ^ + ^jvK- AoVK + \K - l\\ 

n n 

where is the solution of (P) . 



Proposition 1.4 Let us assume that A > and that (HI) to (H4) are satisfied. Let fi be 
a positive Radon measure, and let {us,A^) be a minimizer of J^. Then 



1 A f If 

(1.21) \imsup—Je(ue,As) <- bdfi + - / Vh- AoVh+\h-l 

e^O /iex 2 y 2 J 

Q n 

where h is the solution of 
(1.22) 



— div(«4oV/i) + h = in f2, 
h = 1 on dQ. 



Section II is devoted to the proof of Proposition |]3|. Let (u^, A^) be a sequence of 
minimizers and = curlA^. The energy Ji,{us,Ai,) gives two contributions: inside the 
vortex balls and outside. Thus, first we prove Proposition ITT] where the vortex balls Bi 



with centers Pi are constructed and where the vortex energy is bounded from below. We 
define 

(L23) fis = T—'^diSp^. 



Then, Proposition f?l| implies 



1 /" 1 |2 ^ I log^l 



(1.24) ^ / |VW^>^/i|,. 



"ex ^ I I ' "ex 
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which gives the lower bound inside the vortex balls. The next step is to pass to the limit in 
the energy outside the vortex balls. Letting Hq be the weak limit of hi^/hcx, we obtain 
the following, which is similar to a standard result in homogenization theory 

(1.25) liminf / ^—^> Vho-AoVho. 

J Ciehlx J 

This requires to introduce an auxiliary problem before applying the homogenization theory 
result and it works because the vortex balls are small and thus can be taken out of the 
first integral. 

Finally we derive from the Ginzburg-Landau equations the crucial fact that he satisfies 

(1.26) i- (-div (^) + = + i, 

where ifj^ tends to and /i^ defined in (p..23|) tends to some /io, both convergences being 



strong in W~^''^ for r < 2. The notion of if-convergence and a priori estimates allow us to 
pass to the limit in ( |1.26[ ) in order to get that the weak limit of he/ hex, that we call ho, 
solves 

(1.27) -div{AoVho) + ho = /io- 

Combining the lower bounds of the energy inside and outside the vortex balls ( [1.24| )-( [r^ ), 
we find 

limM ^Je{Ue,Ae) > E{ho) > E{h,). 

The last inequality is true because ( [1.27| ) implies that ho is in V. 

Section III is devoted to the proof of Proposition The proof holds for any positive 
Radon measure /i. We apply it to /i* to get that 

Urn sup -^Je{Ue,Ae) <E{h^), 



e^O hi 



ex 



which will imply the desired results of convergence. 



The upper bound of Proposition |L4| is obtained by constructing test configurations as 
follows. First, given a positive Radon measure /i, we construct approximate measures /i^ 
which converge weakly to /i: 

^ Tie 
t=l 

where /x* is the line element on the circle dB{pl, e) normalized so that /i* (9-B(p* , e)) = 2it. 
The measure fie describes the vortices of our test-configuration. The difficulty is to choose 
the points pi satisfying a number of properties. We tile fl with squares K of size S{e). In 
each square, there is at least a point px where Pe = 0. We choose rij^ points pi regularly 
scattered around px in a ball of radius 1/hex- The number nx is chosen depending on 
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fi{K) so that /ie converge to /i. Once the vortices are constructed, the rest follows easily: 
the magnetic field is defined to be the solution of 



(1.28) 




Then, we are the able to construct a configuration {us,A^) such that curlA^ = and 
has vortices at the points . Moreover, we obtain 

J,K,A,)^1 / l\VK\' + \h,-l\\ 

2 J CLe 

n 

Finally we are able to show that 

limsup / — |V/i,|2 + \he - I? < ^ I h d^i + ]- ! Vh ■ AoVh +\h- Ij^, 
n n n 

where h solves — div(^4oV/i) + h = fi and h = 1 on dQ. 



II Lower bound 

In the following, we will denote V^m = Vm — iAu. We will often drop the subscripts e. 
We consider (m^. A,.) a family of minimizers of Jg, thus a family of solutions of (G.L.). We 
can state a few a priori bounds. Firstly, by the maximum principle, \ue\ < maxa^ < 1. 
Secondly, by minimality, comparing with (ae,0), we get 

Je{us,A,) < J,(a„0). 

But, by hypothesis (H2) on a^, 

Ma„ 0) = l [ I Va,r + 0{hl) < ^ + 0{hl) < Chi. 
n 

Hence, we have the a-priori estimate 

(II.l) Uu,,A,)<Chl. 

In addition, by applying a gauge-transformation to (ms,^^), we can choose the Coulomb- 
gauge divAg = in f2, with A^.n = on dVL. With this choice of gauge, we are easily lead 
(see [^ , |SS1|| ) to the a priori bounds 



(11.2) < CK 

(n.3) ||VMe||L2(Q) < C/le: 

We begin with the proof of Proposition pTT[ 
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II. 1 Proof of Proposition |I.1| 

- Step 1 : Let {u^A) be an energy-minimizer. Denoting \u\ by p, since |Vup > |Vpp, 
we deduce from (|IL1|) : 

(II.4) 



|Vpr + ^(p^-a,)2<CC. 



But, 



|Vp|' = y |V(p->;)P + |Vv/^p-2V(p- v^)- vv^ 
> y |V(p - v^)|2 - 2|V(p - v^)||Vy^|. 
Hence, in view of ( [II. 4| ), 



c 



< C'^ex+^l|V(p-v^)|U2, 



T]{e 



and, since ^ < /^ex, 



1 |V(p - v^)P < m^xiChl, ^) < CC- 



In view of ( p.4| ), we thus have 

(11.5) I J |V(p - + ^(a. - p^)^ < Chi < C\\og e\'. 

n 

- Step 2 : For any t G M, let fit = {x G fi/|p — ^/a^\{x) > t} and 7t = 9fif. Applying the 
coarea formula and arguing as in Lemma IV. 2 of [[SS2||, 



C|Ioger>y"|V(p-v^)r + ^(a,-p2f > |y'|V(p-v^)|K-p2 



> 



C 



r{-it)tdt. 



Here, as in |^S2|] , r(7f) is defined as the infimum over all finite coverings of 7t by balls 
Bi, Bk of the sum ri + ■ ■ ■ + Vk where r j is the radius of Bi . Combining the previous 
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inequality with the mean-value theorem, we find that there exists a t G 
that r(7t) < Ce|log e\^. 



log e| 



such 



- Step 3 : The next step is to construct the vortex-balls : starting from the chosen 7^, 
covered by balls Si, ■ ■ ■ , -B^ (whose sum of the radii is controlled by Cejlog we use 
the method of growing and merging of balls used in Fl, |SS2[| : one needs to grow these 



balls -Bj, keeping a suitable lower bound on the energy they contain, until the desired size 
is reached, with the desired lower bound. When some balls happen to intersect during the 
growth process, they are merged into a larger one. We refer the reader to ||SS2|| , and here 
we only need to apply the result of Proposition IV. 1 of [^S2|] to and f = = e*'^ in 

a = e~v^'^°^ We then obtain the existence of balls Bi = B{pi,ri) such that ( p.l3| ) 
and ( |1.14| ) hold, and 



(II.6) 



I J \Vip-Af + ^J\h- > n\d,\\\og e\{l - 0(1)), 



B,\nt 



B, 



with di = deg(M, dBi) if Bi C Q, and otherwise. But we also have, from the Ginzburg- 
Landau equation — V^/i = p'^i'Vip — A), and from p < 1, 



n n 



p^\Vv-A\^ < / \Vau\^ < Ck 



hence 



\h — hpJ'^ < CrA\h — h. 



exiii4(n) ^ Cri\\h — /iex||/i'i(n) 



Thus, ( |11.6| ) becomes 
(11.7) 



Now, 



< Chl^e~y^\ = o{l). 



A\^>'K\di\\\oge\{l-o{l)). 



1 
2 



Bi\nt 



Bi\Q.t 



> 



> 



Bi\nt 



B,\nt 



mm a. 

B, 



\Vip-AY - 



c 



B^\nt 



|log e\ 



\Vip-A\' 



B,\^t 
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where we have used ( |1.13D . In view of ( |I1.7| ), 

^ j |Vam|' > vr (^mina,^ |rfi||loge|(l-o(l)). 

B^\nt 

So, using the hypotheses (H2) and (H3) on a^, we are led to the two following lower bounds 

(11.8) ^ J IVa^P > nae{pi)\di\\hge\{l-o{l)) 

Bi\nt 

(11.9) i J \\/au\' > 7r6(j9,)|rf,||log£|(l-o(l)). 

Bi\nt 



This proves (|I.15D 



II. 2 Deriving the limiting equation 

For any {pi,di) satisfying ( [1.1 3|) — ([1.15|) , we can define 

27r 

(11.10) J2d^Sp^, 

a measure of vorticity per unit of applied field. We will see that it remains a bounded 
family of measures. 

Lemma II. 1 If A > 0, and {u^, A^) is a family of minimizers of with /i^ = curl A^, we 
can extract a sequence e„ — > such that there exists Hq — 1 & Hq{Q), and fiQ & Ai with 

fi^^ in the sense of measures. 

Proof : As seen in the previous proof, since (m^, A^) is a solution of the second Ginzburg- 
Landau equation 



ex 



and 



n 



Hence, -j^ — 1 is bounded in Hq{Q), and we can find a sequence e„ — > such that ^ 
converges weakly in Hq to some Hq ~ 1. On the other hand, from Proposition [TTI , 



ChJ-^^>J,{u,,A,) > 5^7rM,|6(p,)|log£|(l-o(l)) 



> boY,'^\d^\\^oge\il-oil)), 
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where 60 is given by hypothesis (HI) on Og. Hence, 

2 J ft., - 

n 

thus (/ie„) is a bounded sequence of measures, and extracting again if necessary, we can 
assume that ^e„ converges to some /iq in the sense of measures. • 



Proposition II. 1 Let /iq and ho be the measures and fields defined in Lemma ILL Then 
there exists ro < 2 such that fiQ G W^^'^{Q) Vr G (ro,2), and h^ is the unique solution in 
W^^-" of 



(11.11) 



— div {AoVho) + ho = fiQ in Q 
ho = I on dfl. 



The proof of this proposition requires the following lemma, a slight refinement of the result 
stated in | |SS1| , Lemma II. 3. 



Lemma II. 2 Under the hypotheses of Lemma ILl, for any q > 2, 



hp 



— jje — ^ strongly in {WQ''{Q)y. 



Proof : Denote = i7\ Uj Bi. On Q, \us\ > bo > and = is well-defined. Let q > 2, 
and ^ G Wq''^- We need to show that 



h 



n 



< 



Dropping again some of the subscripts, we have 

(n.i2) 



h 



h 



Then, the method consists in splitting this integral into the integral over the vortex-balls 
(which is going to be negligible because the balls are small enough) and the integral over 
Q, the complement of the balls. 



Step 1 : We prove that 



(n.i3) 



o(l)l|VeiU.(n)- 
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Indeed, since > 60 > 0, 



hex O'e 



< ^^^l|V«|U.(vol(U.B.))i 

Oo "-ex 



where ~ + ~ = ^ and we have used Holder's inequahty twice. Using (|II.3|) 



h 



1 . {^U,Vu) 



< 



In addition, (J^i^f)'' — (Si'^*)'' = '^(1) since we know that ~^ 0- Therefore, ( [II.13| ) 

is proved. 



Step 2 : We observe that 



(11.14) 

We claim that 
(11.15) 

Indeed, 



[iu, V-u) 



1 /■H!(,^^Vt;)-V^e 

''-ex J 0,£ 

n 

[ {zv, Vv) . + 7^ / f — - O i^v, Vv) . V^e 

Aiex J hex J \ as J 



n 



\u\ 



<0{l)\\Vi\\L.. 



\u\ 



- 1 {iv, Wv) ■ W^i 



< 



1 



< c- 



bohcx 
Wv 



{\u\' - a,)\Vv\\V^\ 



lL2(n) 



with 1 + 1 = 1. 

p q 2 



From the a priori estimate (II. 1), 

JiH'- a,Y <C j{\u\'- a,f < Ce'hl = o(l), 
n n 

hence, using || V^'||£^2(Q) < C|| VM||i2(f2) < Chex, we obtain ( [11.151) . Combining ([1I.12|) 
( [11.15| ), we have 



(11.16) 
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Step 3 : We evaluate 



j {iv, Vv) 



V"*"^. Noticing that curl {iv, Vv) = on we have 



an 



dBinn 



There remains to prove that 



(11.17) 



' dBinn 



dv 



in)- 



Let / be a function defined on such that 



bo 



(11.18) 



f{x) = X for X < 
f{x) = 1 for a; > 6o 
|/'(a;)|<C foranyx>0. 



We can define the complex-valued function 
(11.19) w = f{\u\)v. 



It has a meaning everywhere by setting w — u where \u\ < y. Then, it is easy to check 

\Vw\ < C\Vu\ in 



that 
(11.20) 
and 

("-) E/<-|;)-E/<-a 

Using Stokes theorem, we have 
(11.22) 



SB, 



E/v-«.(«».V») + «-«p.))our.(,».V») 



But, on the one hand. 



Bi 



< C&||V«||„|E 



(11.23) 



< o{i)m\\L. 
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as in the proof of ( [il.l3| ). On the other hand, using the fact that, since q > 2, Wq''^ embeds 
in for some /5 < 1, and |curl {iw, Vw)\ < C| Vw p < C| Vwp, we have 



IVmI 



< e 



, ||\77,||2 



(11.24) 



iieii 



o{i)m\ 



since /icx < C|log s\. As in ||SS1|| , the proof remains vahd even if Bi intersects dVt. Combin- 
ing (fT^D, ( PT^ , (|im| ), and ( PT^ , (PTTTD is proved. Consequently, in view of (PTT^ ), 
we can conclude that 



1 



m, Vm) 27r ^ 
h I ~a ^Z.^^^(^^ 



- — / ^ curl ■ 



< 0(1)11^11 



hence that ^curl '■^"^^"^ — ^ strongly in (W^q ''^)' as stated. 

Proof of Proposition II. 1 : For the sake of simplicity, we write e instead of e„. 



- Step 1 : We prove that satisfies 
(11.25) 



he]=fe 



with /e = + ipei where V'e — ^ strongly in iWl''^)' for q > 2. Indeed, we start from the 
second Ginzburg-Landau equation : 



-V^he = {iUe,VA,Ue), 



divide it by and take the curl : 



— div 



curl 



A, 



hence 
(11.26) 



—div 



+ = curl 



+ curl [ AAI 



Now consider a test-function ^ G VFQ''(f2), q > 2, 



ecurl A, 1 



\U\ 



< C||^e||L°°(n)||V^||L2(n)||0£ — |Mn|L2(n). 
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<o(l)l|Ve||L- 



The a-priori bound ( [11. 2| ), || A£||Loo(f^) < 0{hex) and the energy bound, \\a£ — |'uP||l2 < 
Cehex, yield 

ecurl ''''' 

Consequently, curl ^1 — j ~^ strongly in (Wq''^)' for q > 2. Combining this with 
( [11.261 ) and Lemma II. 2, we get the desired result. 

- Step 2 : We prove that fe converges to /xq, the weak limit of /ig, in W~^''^{VL) for any r < 2. 
Indeed, from the upper bound on the energy, we know that ^ V/i^ is bounded in L^(f2), 
hence, in view of ( pi.25[) , is bounded in hence in W~^'^ for p < 2. But, on the 

other hand, fs = fie + i'e, with ip^ bounded in W~^'^ for p < 2, hence fie remains bounded 
in for p < 2. Furthermore, Hs is also bounded in the sense of measures, therefore 

we can apply a theorem of Murat (see [[Mul[| ) which asserts that such a /Xg, bounded in 
the sense of measures and in for p < 2, is necessarily compact in W'^'"^ for r < p. 

Since this is also the case for ipi,, which converges to zero, this implies that is compact 
in W'^'"^ for r < 2. In addition, its limit in the sense of distributions is /io, hence it must 
converge to fio in W~^'^. 



- Step 3 : We wish to pass to the limit in ([II.25[) , but it is not possible directly because the 
i7-convergence requires a right-hand side in H^^. So we are going to pass to the limit in 
the duality sense for a fixed right-hand side. Let g G W~^''^ for q > 2. Using the hypothesis 
(HI) on Qs, (which implies in particular the uniform ellipticity of ^2r), we can apply a 
theorem of Meyers |[Me[| : there exists a. qo > 2, such that if g is in W"^''' with 2 < q < q^, 
then equation 



(11.27) 



-div 



in Vt 

on dVL, 



has a unique solution Ve in W^''^ . Thus, we have 



f 11.28) 



< 



- 1,^ > 



1, and we want to pass to the limit. 



where ^ -|- - 

More precisely, Meyers' theorem yields that the operator which maps g to Ve, is a. 
bounded linear operator from W~^''^ to Wq''^ (for 2 < g < go), hence up to extraction of 
a subsequence, Ve has a weak limit vo in Wq''. We assumed in hypothesis (H4) that j-X 
if-converges to By the definition of if-convergence (see |[MuT[| ), and since Wq''^ C Hq, 
this implies that Vq is the solution of 

—div (AqVvo) + Vq = g in i7 

Vq = on dQ, 



(11.29) 
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Since this possible weak limit vo is unique, the whole sequence Vs converges to vq weakly 
in Wg^''^. In addition, converges strongly to /iq in W~^''^' , thus we have 



w- 



< /e - 1, Ve >wi,g->< /io - 1, ^^0 > • 



On the other hand, — 1 converges weakly to /iq — 1 in Hq. Thus, 



Wr 



l,q' < 



hi 



- >w-^'i—'< ho-l,g > 



Therefore, we can pass to the limit in ([II.28|) , and we are led to 
(11.30) < hQ-l,g >^y-i,g=^_i,,,< fiQ - l,t>o >^i,? . 

Meyers' aforementioned theorem, also yields that for < q' < 2, ( [1I.11| ) has a unique 
solution in W^''^ . Since ( |11.30| ) holds for any g in W~^''^, it implies that is this solution. 



II. 3 Deriving a lower bound outside the vortex balls 

Next, we would like to deduce from ( [11. Ill ) a lower bound like 

liminf [ ^-^^ > f Vho ■ AoVho. 

n\u^B^ n 

But this is impossible to derive straightforwardly because the domain of integration in 
the left-hand side integral is not Q. To remedy this, we replace /i^ by an auxiliary field 
he, a sort of truncated of in the balls. This is a trick that was already used in 
Proposition IV. 1, Step 1. 

Lemma II. 3 There exists h^ such that h^ — 1 E Hq{Q) and 
l)^^-l^ho-ltnHl{Sl), 
2) 

\he-hcy,\ > / h - ftexl ~ o{l) , 



3) 



n\UiB, n 

\Vh 



liminf / ' > / V/iq ■ AV/iq- 



n n 

Proof : We consider a solution of the following minimization problem : 



(11.31) min 

Ae_f/i(n,]R2),div A=o 



j ae\W'^-A\^ + j \cMi\A-h 



2 

ex I ; 
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where V^p denotes the gradient of the phase of Ue which is well-defined in f2\ Uj Bi. If we 
write he = curlA^, and we test ([1I.31|) with hs, we have 



(11.32) 



e\ ' I I ""ex 



In addition, and satisfy the following equations : 
(11.33) 



-V-L/i^ = a^iVif - Ae) in n\ Uj Bi 
he = est = Ci on Bi,^i 

on d^l. 



'^e^ex / hgx 



he ^ex 

Thus, it satisfies 
(11.34) 

where Ue is the measure defined by 

(11.35) Ve e <'^(fi), (g > 2), / = E /T / + 5Z /T / 

On the other hand, using Cauchy-Schwartz inequality. 



B, 



^ ^ he^ 



L2 



In view of (|L32D 
Hence, 



L2 



is bounded, and (^jT'j)^ < Xli'^* ~* from Proposition |LT 



B, 



o{im\\ 



L°° ■ 



On the other hand, the same proof as for Lemma II. 2 shows that 



dBi 



o(l)||ell 



Hence, in view of ( [II.35|) , Ue — fie converges strongly to in {Wq''^)'. The same argument 
as in Proposition II. 1 allows to conclude from ( |11.34| ) that 
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using the uniqueness of the solution of ( |11.11| ). 
Using (IlId and (|IL33| ), we get 



+ \he - he 



n\UiBi 



n\\JiBi 



As in the proof of Proposition we have 



+ o(l). 



Q.\JiBi 



n\UiBi 



Thus, assertion 2) is proved. In addition, — 1 is bounded in Hq^H) and the convergence 
to /lo — 1 is weak in Hq. There remains to prove the third assertion. But it is a classical 
result in homogenization theory (see ||J KO| ] ) that, since — 1 ^ /iq — 1 in Hq^VI) and 
if-converges to A^, 

2 



lim inf 

£-♦0 



V 



hr 



This completes the proof of the lemma. 

We recall that we defined E in ( p9|) . 
Lemma II. 4 With the same notations, 



lim inf 



^l/Uol + ]^jyho- A^Who + |/io - ir = ^(/io). 



Proof : The energy can easily be bounded from below as follows, splitting between the 
contribution inside the vortex-balls and the contribution outside : 

JeiUsAe) > \ j \yAu\^ +\h-h,^\^ 

n 

> \ j \^AU? + \ j p'\V^-A\^ + \j\h-K^\\ 

UieiBi n\u,Bi n 

As previously, since for the energy-minimizers — V^/i = (iw, V^m), and — a^] < 
in Q\ Uj Bi, we have 

|V/i|2 



-(1-0(1)). 



n\UiBi 



n\UiBi 
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Therefore, in view of Proposition [0] , 

Uu„ A,)>nJ2 \d^\b{p^)\log e\{l - o(l)) + f 
i ^ 

and with assertion 2) of Lemma II. 3, 



\Vh\ 



:i-o(i 



JeiUe^A) ^ 1 |log e\ 



2 /i. 



+ 



n n n 

We thus obtain, using assertion 3) of Lemma II. 3 that 



oil] 



(11.36) hminf > liminf- ( / 



hL 



2 /i. 



V/io ■ AoVho + \ho - 11 



Similarly, using ( [II. 8| ), we obtain 
(11.37) liminf:^ife^>liminfl^'^°^" 



2 I K. 



n / n 



Then, using the weak convergence of fis to /iq in A^, and the weak lower semi-continuity 
of /i 1-^ b\fi\, we conclude from ( [1I.36| ) that 



hi. 



The final convergence result will then follow from the combination of this result with 
the upper bound of Section III, leading to the fact that necessarily ho has to be h^, the 
minimizer of E, and /iq = /i*. 

Ill Upper Bound 

In this section we prove Proposition 1.4. First we remark that if h is the solution of 
—div{AVh) + h = n with boundary value 1, then 



h{x)~l= / G{x,y)d{^-l){y), 
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where G{.,y) is the solution of —div{AVh) + h = 6y vanishing on dfl and /i — 1 denotes 
the difference between the measure /i and the Lebesgue measure in Q. From this it follows 
easily that 

(III.l) Jvh-AVh+\h-l\^ = Jj G{x,y)d{fx-l){x)d{fx-l){y). 

n 

This last expression will be the one we use. 

To prove Proposition 1.4 we will then need some properties of the Green functions Gs, 
Go associated to the operators — div(^eVM) +u and — div(.4oV-u) +u respectively. These 
properties will be proved at the end of this section. 

Lemma III.l Let = h + 13;, he a sequence of functions satisfying (HI) to (H4), and 
.Ao be the homogenized limit of the matrices A-e = a^^^X as e goes to zero. For any 
y E Q, let Gs{.,y) (resp. Go{.,y)) be the solution o/ — div(^eVG'e) + Gs = 6y (resp. 
— div(^VGo) + Gq = 6y) that vanishes on dQ. 
The following properties hold: 

1 ) Gsix, y), Gq{x, y) are positive functions, and symmetric in x and y. 

2) A denoting the diagonal in R^, there exists C > such that Ge{x,y), Go{x,y) are 
bounded by 

C{\\og\x-y\\ + l) 

for all x,y e Q X Q \ A. 

3) For any compact K dVt, there exists C > such that for any x,y E Q 

Ge[x, y) + — — log \x-y\< 



27r ^'-ri{e)' 

where rj{e) is defined in (H3). 

4) Gs converges to Gq locally uniformly in Q x Q\ A. 

Then we have the following easy Lemma: 
Lemma III. 2 The function 

(in.2) ^(/^) = ^ jbdfi+^ jj Go{x,y)d{fi-l){x)d{fi-l){y) 

is sequentially lower semicontinuous over the set of positive Radon measures supported in 
fl, with respect to weak-* convergence. 



The proof of this can be found in |^] for instance. Note that /(.) is well defined over the 



set of positive Radon measures if we admit the value +oo. Note also that if we restrict 



to measures in H-'^{n) then (|ITLTD shows that /(/i) is a lower semicontinous functional of 
h = L^^jj, where L is the operator u — div(^V-u) + u defined on HliVL). It follows that 
J is a lower semicontinuous function of /i with respect to convergence. 
Now the proof of Proposition 1.4 splits into two propositions. 



26 



Proposition III.l Assume that A > and that (HI) to (H4) are satisfied. Let fi be a 
positive Radon measure with support in Q and {pl)i<i<n^ be families of points in Q such 
that Vi 7^ j 

(IIL3) \pl -p^^\> 4e, d{pl, dQ) >ao>0, 

where is independent of e, 



(III.4) 

and 

(III.5) 



fi, in the sense of measures 



i=l 



( 



lim 

e^O 



E 



\|pi:-pil<" 



I log be 

hi. 



J 



Then there exist configurations {vi;,B^)^yQ such that 

riTT a\ r Mve,B,) A 

[iii.Dj limsup — < — limsup 

£^0 hf. 2 /lex 



0. 



Go d{fJ' — l)d{fi — 1), 



where Gq is defined in Lemma IILl. 



This proposition states that under reasonable hypotheses on points , one can construct 
a good test configuration with prescribed vortices at p\. Moreover, ( [111.4] ) implies that 
ne/hex is bounded. The following Proposition asserts that the construction of points is 
possible. 

Proposition III. 2 Assume that A > and that (HI) to (H4) are satisfied. Then given 
any positive Radon measure /i of the form a{x) dx where a is a positive continuous function 
compactly supported in fl, there exist families of points {pl)i<i<n^ satisfying ( [1II.3|) , (|III.4|) , 
( [1II.5I) and such that 

(III. 7) limsup < / b{x)dn{x). 

£^0 hex J 

n 

The proof of Proposition 1.4 follows easily from these two Propositions. First, taking 
any positive Radon measure fi supported in fl, we may approach it in the weak-* topology 
by measures /i„ = o"„(a;) dx where cr„ G Cc{fl) is a positive function. Applying Propositions 
III.l and III. 2, we may construct test-configurations (f",-B")e>o such that 
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J (v^ B^) A. r 1 

limsup — — < — b{x) + 7: I I Go d{^n - l)c?(/in - 1) 



ex 



hi. - 2 ' ' ' 2 



Therefore the same inequahty is satisfied if we replace (f -B") by the minimizing config- 
uration {ue, As). This proves that for each n, 



hmsup J^' < /(/in), 



and then, using Lemma III. 2, 



(IIL8) hinsup ^^^^^^ < ^ Jbdf,+ ^ 1 1 Go{x,y)d{^,-l){x)d{^,-l){y). 



n 



Using ( |III.1|) we get the conclusion of Proposition 1.4. 
III.l Proof of Proposition III.l 

The method for constructing a test configuration [vs, -Be) with prescribed vortices {pl)i<i<n^ 
follows closely that of [[SS3||. First we define to be the solution of 



(III.9) 



-div{AeVhe) + he = ^ni IB. Q 

i=l 

hr = /lex on dfl, 



where nl is the line element on the circle dB{pl, e) normalized so that fil{dB{pl, e)) = 27t. 

Then we let B^ be any vector field such that curl B^, = h^. Finally, we define = PsG^'^" 
as follows: first we let 



(III.IO) pe{x) 



if \x — \ < e for some i, 

a/ ae(x)- — — - if e < |x — I < 2e for some i, 

\/aJx) otherwise, 
and for any x ^ Qe = ^ \ ^iB{pf, e), 



(III.ll) <f,{x) = j (5, 



A,V^K).Td^, 



where Xq is a base point in $7^, (xo,a;) is any curve joining xq to x in fi^ and r is the 
tangent vector to the curve. From (pi.9| ), we see that this definition of ^Pe^x) does not 
depend modulo 2tt on the particular curve {xq^x) chosen. The fact that <^e is not defined 
on UiB{j)l,e) is not important since pe is zero there. Thus, ips satisfies 

(III.12) -A,V^K = Vp^e-B, 
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in Qs- Having defined Vs = PeG^'^^ we estimate Je(fe, 5e). Recall that 

(111.13) J,{V„B,) = ^ j |Vper + Pe'|V(/^e-S,|2+|/i,-/i,,|2 + -L(a,-(^,2^)2_ 



Using the fact that | Vae| <^ /lex (hypothesis (H2)) and that the number of points p\ is less 
than Chex — which follows from ( [1II.4|) — it is not difficult to check that 



(III.14) 



2 

ex' 



f2 



Also, from ( prToi) , (|ITT:T2D , 

in fi^. Therefore, replacing in ( |111.13| ) and in view of ( |111.14| ) 



(111.15) 



hm sup — < hm sup —— 



Because is the solution of ( [111.9| ), we may rewrite the right-hand side of this inequality 



as 



where 
(III.16) 



limsup- // Geix,y)d{p^-l){x)d{pe- l)iy), 



and p\ is defined in ( [TTLQD . It follows from ( [TTLl ), (|Tr|) and ( ITTLTel) that /i^ ^ p as 
e ^ 0. Thus, to finish the proof of the proposition, it remains to show that 



(III.17) limsup - // Ge (i(p£ — 1) (i(/i£ — 1) < ^ limsup 

e-+0 



27r J]a,(p^ 



4 = 1 



hp 



+ - Go d{p-l)d{p-l] 
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Proof of ( pmD 

Let a > and let Aq = {{x,y) \ \x — y\ < a}. Recall that fi^ —>■ Hence, it follows that 
(yUe — 1) ® (/ie — 1) ^ (/i — 1) (8> (/^ — 1) as £ ^ 0. But from Lemma II. 1, tends to Go 
uniformly in x i7 \ A^, therefore 



(III.18) 



lim - 

£^0 2 



Ged{fls - 1) d{fle - 1) 



1 



God{fi-l)d{fi~l). 



Now we treat the integral on Aq,. More precisely we prove that 



rie 



(111.19) lim sup / / G, d{fie - 1) d{fie - 1) < - lim sup 

2 e^O hex 



where lim„^oOa(l) = 0. Adding (|TlT|), (|TlT|) and letting a^O yields ( pUTl) . We are 
left with proving ( [1II.19|) . First we use the bound \Gs{x,y)\ < C |log |x — ?/|| from which 
one easily gets 

J J Ge d{fle — I) d{He — 1) < J J (i/ie (i/i^ + Ctt^ | log « | . 
Ac Ac 

Therefore ( [1II.19| ) will follow if we prove 



Tie 



(III.20) 



lim sup 1 1 Gs dus dfi^ < — lim sup 



£^0 



271 J2(^s{pl 



i=l 



hp 



Ac 



To prove this, we come back to the definition of /i^. From this definition, we have 



(III.21) 



Ge dlle dfle < 



hi 



\\pl-pi\<2a 



i=l 



I 



Let us first estimate the first sum on the right-hand side. If a; G Supp \i\ = dB{pl,e), 
y G Supp fii and i ^ j, since |p* — pH > ie, then \x — y\ > — Using the bound 
\Gs{x, y)\ < G |log |x — I together with the fact that \pl —pl\ < 2a and a is small enough, 
we get 

rf .... 
Gedfildfil < G log|p^ -p^l 
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Then, by hypothesis ( |111.5| ) 
fIII.22) 



hmsup-^ 

\pi-pi\<2a 



Gs dfil dfil < Oq,(1). 



As for the second sum in the right-hand side of ( [III.21|) , we use property 3) in Lemma III.l 
to get that for any I < i < Us, and any x,y E Supp /x* , 



(III.23) 



Ge[x^y) H — log \x — y\ < ^ |log e\ 



27r ~ ri{e) 

But X G Supp/i* is equivalent to \x — p\\ = e. Then property (H2) of implies that 



a^{x) ~ a^{p\) as e — > 0. Replacing in ([III.23D and integrating w.r.t. yU* ® nl yields 



Ge di^l d^i\ < 27raM) |log ^1 (1 + o.(l)) 
and then, summing over 1 < i < and dividing by /lex, 

Tie 

27r^a,(p:) 



(111.24) 



limsup — ^ 



.0 hj^ 



1=1 



Gs d^l dfil < — lim sup 

2 e^O 



i=l 



Here we have used the fact that |log e\ ~ Ah^x- Thus ( [III.2CI| ) is proved and the Proposition 
follows. • 



III.2 Proof of Proposition III.2 

Let /i = (7{x) dx, G = \\u\\oo and ao = dist{suppfi,dQ). Also, let 

(IIL25) n = {xen\ d{x, on) > ao/2}. 

Recall that from hypothesis (H3) on there exists a positive function 6{6) such that 

1 



(IIL26) 



(log |log e|) 



1 ; 

2 



and for any x E fl, min 6^ = 0. 

B{x,S{e)) 



For any e > 0, we tile with open squares of sidelength 26 (e) and let IC{e) be 
the family of those squares that are entirely inside fl. We denote by ck the center of a 
square K. Since /i is absolutely continuous with respect to the Lebesgue measure, we have 
fx{K) < G5\ 

Now the family of points {pl)i<i<n^ is defined as follows: for any K G /C(e), we let 



(IIL27) 



n{K, e) 



2ti 
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where [x] is the biggest integer no greater than x. Using ( |ill.26| ) there is a point pk G 
B{ck,S) such that P^ipx) = {px is a pinning site). We now pick n{K,e) points evenly 
scattered in the ball B{pK,l/hex), and we call V{K,e) their union. By evenly scattered 
we mean that for any p,q E V{K, e), 



(III.28) 

We let 
(III.29) 



\p-q\> 



c 



hexVn{K,e) 



KeK{e) 



and Vie) = UKeK{e)V{K, e) = {p\ 



be our family of points. We now check that this family satisfies ([111. 3D , ( [111.4] ), ( [111.5[) and 

(io). 

( [111.3| ) is clear from ( [111.28[ ) if belong to the same pinning site. It is even more 
true if p\,p^ir do not belong to the same site since in this case their mutual distance is at 
least 28{e) > e. Moreover from ( [111.25[ ) we have d{pl, dVt) > ao/2. 

For (pL4[) , let 



(III.30) 



27r 



i=l 



and / be a continuous function in Q. We let 7^ = sup;^£/c(e) snp^y^^^ Ifi^) — fiy)\- Then 
since the size of the squares in JC{e) tends to zero with e, so does 7e. Let be the union 
of the squares in /C(e), then for e small enough supp/i C and 



< 



^ J2 |M^)-/X.(i^)|+7e(/Xe + /x)(i^e)- 
Ke!C{e) 



It is clear that the second term on the right-hand side goes to zero with e. For the first 
term we note that from (PI.27[) , ([III.30[) , we have |/i(-ft") — fis{K)\ < 27r//iox while the 
number of squares in /C(e) is of the order of 1/6"^. From ( [III.26[) it then follows that 
J2k&jc{s) f^e{K)\ tends to zero with e. We thus have lime_^o / f dfie = J f dfi and 

follows. 

We easily deduce ( [III.7[ ) from ( [III.4[ ). Indeed from (112) and the fact that each point 
is at a distance at most l/h^x from a pinning site, we get that a^^p) ~ b{p) as £ — ^ 0, 
uniformly in p EV{e). Moreover, since rie/hex is bounded. 



2nJ2aeipi 



2n 



lim ■ 

£^0 



i=l 



hp 



lim ■ 



i=l 



h{p\ 



he. 



b{x) dfi{x), 



by the convergence of /i^ to /x. 
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It remains to prove ( p.ll.5| ). We split the sum in ( [I11.5| ) as follows: let T{e) be the set 
of pairs of indices (i, j) such that 1 < i 7^ j < and p*,^^ belong to the same square of 
the subdivision lC{e). Let J{e) be pairs such that belong to different squares. 

Then 

(111.31) |iogb:-pii|= |iogb:-pii|+ |iogbe-i4i| 

\pl-p'e\<a \p\-pi\<a \pl-pi\<a 

The first sum in ( |111.31| ) is estimated as follows. For every K G lC{e), fi{K) < C6^ thus 
the number of points of V{e) in K is less than CS^hex- The number of squares being of the 
order of 6-^, the cardinal of X(£) is less than CS^hl^. Using (pi^ , { ^T?27\ ) and ( PI^ ), 
we find 

(111.32) Y |logb^ -Pill < ChlSHog |log £| « hi. 
(ij)ex(£) 

To treat the second sum in ( pil.31|) , we note that if K and K' are distinct squares in 
}C{6) and p & K, q & K' then 

\/x e K,\/y e K', \x - y\ < A\p - q\. 

Thus we may write, using the fact that fJ.{K) < Cd"^ 



Y |logbe - pill < Chl^ {\\og\x-y\\ + l)dxdy. 



^"^i KxK' 
pl€K,pieK' 

Summing over pairs of squares K, K' G IC{e) such that KxK' intersects {(x, y) \ \x — y\ < 
a} we get for e small enough 

(III.33) Y H\pl-p'e\\<Chl Jj {\\og\x-y\\ + l)dxdy. 

(i'j)&J{e) \x-y\<2a 
\pl-pi\<a 

Summing (|III.32|) , (pi.33| ), dividing by h^^ and letting e and then a tend to zero yields 
( [nrsD . Proposit ion III. 2 is proved. • 



III. 3 Proof of Lemma III.l 

The fact that and Gq are positive is a simple consequence of the maximum principle, 
that they are symmetric is standard and follows from Green's identity. 
The inequality 

Geix, y), Go{x, y) < -C log \x-y\ + C 



33 



is a well known property of Green functions for elliptic operators in divergence form, a 
proof can be found in [pt| . 

To prove property 3), we let 

v,{x, y) = Geix, y) + log \x - y\ 
and be the operator u i— > — div(^eVM) + u. Then letting = Li,Vi;{-,y), we have 

(III.34) y) = -^V^. V. log \x-y\-^ log |x - y|. 

Thus for any 1 < q < 2, there is a C independent of y and e, such that 2/)||l9 < 

C/ri{e). On the other hand, Vs{-,y) is bounded in W^''^{Q) independently of e and y (see 



Now, Theorem 2 of |[Me|| implies that there exist p > 2 and p' < 2 such that if m satisfies 
L^u = /, then for any compact K G Q, 



We may choose q < 2 such that W ^'^ C L'^ and p' < q. Thus, we find that v^{.,y) is 
bounded in W^'''^{K) by C/ri{e). Since p > 2, this yields the uniform bound Vx E K,\ly E fi, 

r]{e) 

i.e. property 3). 

To prove property 4), we note that for any a > 0, L^Gsi-, y) = in Q \ B{y, a) while 
Gs{-,y) is bounded in W^'''{Q) independently of e and y (see [^). Using the aforementioned 
result of ipviell , we find that G^{.,y) is bounded in Wl^^iVt \ B{y,a)), for some p > 2, 
independently of y and e, thus G^ converges locally uniformly in O x f2 \ A, where A is the 
diagonal. The limit is necessarily Gq, since Go{.,y) satisfies LQGo{.,y) = — div^VajGo + 
Go = and i^-converges to Lq. Lemma II. 1 is proved. • 

IV Convergence results 

We can then proceed as in the rest of Section III in [|SS3 



Proposition IV.l The minimum of E is uniquely achieved by G C^''^{Q)(i'~f < 1] 
satisfying 



(IV.l) 



> 1 in il 

h^: = 1 on dQ 

/i* := -div{AoVK) + K>0 

K- I1-— ll/i* = 
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As in ||SS3|| , we divide the proof of this proposition into several lemmas. 



Lemma IV. 1 Let fi^ and fi^ be the positive and negative parts of the measure /i*. Then 

/i* = 1 - — (xj a.e. 
Ab 

Ai* = 1 + — yu^ a.e. 

Ab , Ab 

1 <K<1 + — . 

2 - 2 

Proof : As in |[SS3|| , the minimum of E is achieved by some h.^, by lower semi-continuity. 
Performing variations (1 + tf)fi^, where / G C^{Q), and looking at the first order in t ^ 0, 
we find similarly as in [ [SS3| that 



2 

Hence 



Ab, , 

— + [K — = 0. 



Ab 

K = l- — nJ a.e. 
ti^ = 1 + — a.e.. 



As in [ ^S3|| , considering variations /x^, + u, where u E Ai Ci H ^ and u and /i* are mutually 
singular, we are led tol — ^</i*<l-|-^. • 

Lemma IV. 2 yU^, is a positive measure. 

Proof : 



fi^{h^ - l)^ = j fi^{h^ - 1)+ - y fi^ - !)+• 
Since — 1)+ = /i^-a.e., we have 

j - 1)+ = - j - 1)+ 

n n 

= j VK-{AoVh,) + h,{K~l)>0, 

h,>l 
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because ,Ao is a symmetric positive matrix (this follows from the compactness of the set 
of matrices bounded from above and below). We deduce that 

- 1)+ = 0, 



n 



but since /i* — 1 = /i^ a.e., we have 

Ab 



h 

hence /i~ = 0, and yU* > 0. • 

Thus, /i* satisfies all the properties listed in ( [IV. 1|) . 

We can now complete the convergence results. From the upper bound of Proposition 
1.4 and Lemma II. 4, we deduce that for our family of minimizers {us, A^), 

minE = EiK) > liminf > E(ho) > EiK). 

V ''^^ 

being the unique minimizer of E, we conclude that = and thus /io = /i*- We also 
obtain 

(IV.2) lim-%^ = i?(M. 

^ ex 

Since the possible limits are unique, the whole family converges to h^, and the same 

for fXe- 

In view of ( p.37|) , we have 
lim inf ^ — — > lim inf - ^ 



^0 /iL <^^0 2 \ /lex 



Q / n 



while 



> ^ y b\f^*\ + \ j V/i* ■ A^WK + \K - l| 



limsup < 4 / ^^^*\ + \ I ^K- AoVK + \K - ip 



ex 



Thus, we deduce that 



lim J a^l/iel = j bfx^. 
n n 
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On the other hand, 



hminf / ae|/ie| > hminf / > / 
J '^^o J J 



hence b\l^e\ Jq bfi^, while bfie bfi^. We conclude that b{\fis\ — fis) and 
thus l/igl and fis have the same limiting measure This proves ( |I.16D , (tI.17|) , and ([I.18| ). 
Following |[SS3|| , Section IV, we can also prove easily the following : 



Proposition IV. 2 If A = 0, then K = I and ^-1^0 strongly in HI{VL). If A > 0, 
then — 1 ^ — 1 in Hq{Q), the convergence is not strong and 

^—^ > VK ■ AqVK + A&/i* in M. 

Proof : First, it is easy to get, as seen in Lemma II. 4 for example, that 



|Va,m.|' > 



(1 -o(l)). 



lim inf 



hl 



thus, we have 
(IV.3) 
(IV.4) 



The case A = follows easily from the upper bound min Jeiug, A^) < o^h"^^) of Section II 
combined with ( [IV.4|) . 

The convergence of ^ to /i* is weak in H^, in general, thus strong in L^(f2), and 




lim 



he 

hr 



- 1 



\K - iP. 



ex 

n n 



Combining this to the convergence result ( [IV. 2|) , we have 
(IV.5) 



lim - / 

6^0 2 J hta, 
n 



2 a/" 1 



n 



n 
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Then, we argue as in ||SS3|| , Proposition IV. 1. Roughly speaking, one considers any open 
set U C Q, and gets a lower bound 

limmt / — = limmt / — h / 

u un{u,B,) u\u^B, 



> A y b\l2e\ + J VK- AoVK 
u u 

> A J bix^ + J VK ■ AoVK. 
u u 

Since this is true for any f/ C fi, comparing this to ( |1V.4| ) and ( |1V.5| ), we obtain as in 

|2 



|SS3 



VK ■ AoVK + Abjj.^ in M. 



This completes the proof of Theorems 1, 2 and 3. 
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